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1 Introduction 



Measurement of mixing-induced CP violation in B° s decays is of prime importance in 
probing physics beyond the Standard Model. Final states that are CP eigenstates with 
large rates and high detection efficiencies are very useful for such studies. The — > 
J/ip fo(980), /o(980) — > 7r + 7r~ decay mode, a CP-odd eigenstate, was discovered by the 
LHCb collaboration [l] and subsequently confirmed by several experiments [2j. As we use 
the J/ip ^ + fi~ decay the final state has four charged tracks, and has high detection 
efficiency. LHCb has used this mode to measure the CP violating phase S [3], which 
complements measurements in the J/ip(j) final state |4j[5]. It is possible that a larger 7r + 7r~ 
mass range could also be used for such studies. Therefore, to fully exploit the J/ipn + it~ 
final state for measuring CP violation, it is important to determine its resonant and CP 
content. The tree- level Feynman diagram for the process is shown in Fig. [T| 




Figure 1: Leading order diagram for B° s decays into J/ipn + n . 

In this paper the J/ifm + and 7r + 7r~ mass spectra, and decay angular distributions 
are used to study the resonant and non-resonant structures. This differs from a classical 
"Dalitz plot" analysis |6] because one of the particles in the final state, the J/ip, has 
spin-1 and its three decay amplitudes must be considered. We first show that there are 
no evident structures in the J/ipn + invariant mass, and then model the 7r + 7r~ invariant 
mass with a series of resonant and non-resonant amplitudes. The data are then fitted 
with the coherent sum of these amplitudes. We report on the resonant structure and the 
CP content of the final state. 



2 Data sample and analysis requirements 

The data sample contains 1.0 fb" 1 of integrated luminosity collected with the LHCb detec- 
tor (7| using pp collisions at a center-of-mass energy of 7 TeV. The detector is a single-arm 
forward spectrometer covering the pseudorapidity range 2 < 77 < 5, designed for the study 
of particles containing b or c quarks. Components include a high precision tracking sys- 
tem consisting of a silicon-strip vertex detector surrounding the pp interaction region, a 
large-area silicon-strip detector located upstream of a dipole magnet with a bending power 
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of about 4 Tm, and three stations of silicon-strip detectors and straw drift-tubes placed 
downstream. The combined tracking system has a momentum resolution Ap/p that varies 
from 0.4% at 5GeV to 0.6% at 100 GeV (we work in units where c = 1), and an impact 
parameter resolution of 20 [xm for tracks with large transverse momentum with respect 
to the proton beam direction. Charged hadrons are identified using two ring-imaging 
Cherenkov (RICH) detectors. Photon, electron and hadron candidates are identified by 
a calorimeter system consisting of scintillating-pad and pre-shower detectors, an electro- 
magnetic calorimeter and a hadronic calorimeter. Muons are identified by a muon system 
composed of alternating layers of iron and multiwire proportional chambers. The trig- 
ger consists of a hardware stage, based on information from the calorimeter and muon 
systems, followed by a software stage which applies a full event reconstruction. 

Events selected for this analysis are triggered by a J/ip ^ decay. Muon candi- 

dates are selected at the hardware level using their penetration through iron and detection 
in a series of tracking chambers. They are also required in the software level to be con- 
sistent with coming from the decay of a B® meson into a J/ip. Only J/ip decays that are 
triggered on are used. 

3 Selection requirements 

The selection requirements discussed here are imposed to isolate B° s candidates with high 
signal yield and minimum background. This is accomplished by first selecting candidate 
J/ip — > fi + fi~ decays, selecting a pair of pion candidates of opposite charge, and then 
testing if all four tracks form a common decay vertex. To be considered a J/ip — > fi + fi~ 
candidate particles of opposite charge are required to have transverse momentum, p^, 
greater than 500 MeV, be identified as muons, and form a vertex with fit x 2 P er number 
of degrees of freedom (ndf) less than 11. After applying these requirements, there is a 
large J/ip signal over a small background |1|. Only candidates with dimuon invariant 
mass between —48 MeV to +43 MeV relative to the observed J/ip mass peak are selected. 
The requirement is asymmetric because of final state electromagnetic radiation. The two 
muons subsequently are kinematically constrained to the known J/ip mass (8j. 

Pion and kaon candidates are positively identified using the RICH system. Cherenkov 
photons are matched to charged tracks, the emission angles of the photons compared with 
those expected if the particle is an electron, pion, kaon or proton, and a likelihood is then 
computed. The particle identification is done by using the logarithm of the likelihood ratio 
comparing two particle hypotheses (DLL). For pion selection we require DLL(7r — K) > 
-10. 

Candidate 7r + 7r~ combinations are selected if each particle is inconsistent with having 
been produced at the primary vertex. This is done by use of the impact parameter (IP) 
defined as the minimum distance of approach of the track with respect to the primary 
vertex. We require that the \ 2 formed by using the hypothesis that the IP is zero be 
greater than 9 for each track. Furthermore, each pion candidate must have px > 250 
MeV and the scalar sum of the two pion candidate momentum, Pt{^ + ) + Pt^'), must 
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m(JApK + K') (MeV) m(K + K") (GeV) 



Figure 2: (a) Invariant mass spectrum of J/ipK + K~ for candidates with m(K + K~) < 
1050 MeV. The data has been fitted with a double-Gaussian signal and linear background 
functions shown as a dashed line. The solid curve shows the sum. (b) Background 
subtracted invariant mass spectrum of K + K~ for events with m(K + K~) < 1050 MeV. 
The dashed line (barely visible along the a:-axis) shows the S-wave contribution and the 
solid curve is the sum of the S-wave and a P-wave Breit-Wigner functions, fitted to the 
data. 



be greater than 900 MeV. To select B° s candidates we further require that the two pion 
candidates form a vertex with a \ 2 < 10, that they form a candidate B° s vertex with the 
J/ip where the vertex fit xV n df < 5, that this vertex is greater than 1.5 mm from the 
primary vertex and the angle between the B Q a momentum vector and the vector from the 
primary vertex to the B° s vertex must be less than 11.8 mrad 

We use the decay B° s — > J/ip<f), <p K + K~ as a normalization and control channel 
in this paper. The selection criteria are identical to the ones used for J/ifjir + 7T~ except 
for the particle identification requirement. Kaon candidates are selected requiring that 
DLL(K-n) > 0. Figure[2|a) shows the J/tpK + K~ mass for all events with m(K + K ) < 
1050 MeV. The K + K~ combination is not, however, pure due to the presence of an 
S-wave contribution |9] . We determine the yield by fitting the data to a relativistic P- 
wave Breit-Wigner function that is convolved with a Gaussian function to account for the 
experimental mass resolution and a straight line for the S-wave. We use the ^PZot method 



to subtract the background 10 . This involves fitting the J/ijjK + K~ mass spectrum, 
determining the signal and background weights and then plotting the resulting weighted 
mass spectrum, shown in Fig. |2](b). There is a large peak at the <fi meson mass with a 
small S-wave component. The mass fit gives 20,934±150 events of which (95.5 ± 0.3)% 
are <fi and the remainder is the S-wave contribution. 

The invariant mass of the selected J/ifm + 7r~ combinations, where the dimuon candi- 
date pair is constrained to have the J/ip mass, is shown in Fig. [3] There is a large peak at 
the B° s mass and a smaller one at the B° mass on top of a background. A double-Gaussian 
function is used to fit the signal, the core Gaussian mean and width are allowed to vary, 
and the fraction and width ratio for the second Gaussian are fixed to that obtained in the 
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Figure 3: Invariant mass of J/iIjtt + tt^ candidate combinations. The data have been fitted 
with double-Gaussian signal and several background functions. The (red) solid line shows 
the B° s signal, the (brown) dotted line shows the combinatorial background, the (green) 
short-dashed shows the B~ background, the (purple) dot-dashed is B° — > J/tpn+n^, the 
(black) dot-long dashed is the sum of B° s — > J/iprj' and B° s — > J/ip4> when — > 7r + 7r~7r° 
backgrounds, the (light blue) long-dashed is the B° — > J/i/jK~7t + reflection, and the 
(blue) solid line is the total. 



fit of B° s — > J I ipcf). Other components in the fit model take into account contributions from 
B~ ->• J/ipK-(n-), B° s ->• J/HiV' -> PI, B° s J/if><j>,</> tt+tt-tt , B° J/ipn+iT- 
backgrounds and a B° — > J/if)K~ir + reflection. Here and elsewhere charged conjugated 
modes are used when appropriate. The shape of the B° — > J/ipir + TT~ signal is taken to 
be the same as that of the B®. The exponential combinatorial background shape is taken 
from wrong-sign combinations, that are the sum of 7r + 7r + and 7t~tt~ candidates. The 
shapes of the other components are taken from the Monte Carlo simulation with their 
normalizations allowed to vary (see Sect. 4.2). The mass fit gives 7598 ± 120 signal and 
5825 ± 54 background candidates within ±20 MeV of the B° s mass peak. 



4 Analysis formalism 

The decay of B° s — > J/ilm + ix~ with the J/ip — > fi + ^~ can be described by four variables. 
These are taken to be the invariant mass squared of J/ip7t + (s^ = m 2 (J/ijjn + )), the 
invariant mass squared of tt + tt~ (S23 = m 2 (7r + 7r~)), the J/ip helicity angle (6j/^), which is 
the angle of the fi + in the J/ ip rest frame with respect to the J/ip direction in the B° s rest 
frame, and the angle between the J/ijj and 7r + 7r~ decay planes (x) in the B° s rest frame. 
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Figure 4: Background subtracted x distribution from B° s — > J/ip-n + -n~ candidates. 

To improve the resolution of these variables we perform a kinematic fit constraining the 
B° s and J/tfj masses to their PDG mass values (8), and recompute the final state momenta. 
To simplify the probability density function (PDF), we analyze the decay process after 
integrating over x, that eliminates several interference terms. The x distribution is shown 
in Fig. [4] after background subtraction using wrong-sign events. The distribution has little 
structure, and thus the x acceptance can be integrated over without biasing the other 
variables. 



4.1 The decay model for B® — >• J/?/?7r + 7r 

One of the main challenges in performing a Dalitz plot angular analysis is to construct 
a realistic probability density function (PDF), where both the kinematic and dynamical 
properties are modeled accurately. The overall PDF given by the sum of signal, S, and 
background, B, functions is 

F{si 2 , S 23 , 0.]/^) = ^^S(S12, S 23 , 9j/^)S(s 12 , S 23 , 6j/i>) + — ./ S1 ^ B(S 12 , S 2 3, 9 J/J,), (1) 

yVsig yVbkg 

where / S j g is the fraction of the signal in the fitted region and e is the detection efficiency 
The normalization factors are given by 



•A/rig = J £{s 12 ,S2 3 , 9j/tp)S(s 12 , s 23 , 9 j/^)ds 12 ds 2 3d cos 9 J/in 
A/"bk g = J B(s 1 2,s 2 3,0 J /^) ds 1 2ds 2 3dcos9j/. l p. (2) 
In this analysis we apply a formalism similar to that used in Belle's analysis of B° — > 



K tt + Xci decays 11 
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To investigate if there are visible exotic structures in the J/ipir + system as claimed 
in similar decays 12 , we examine the J/tpir + mass distribution shown in Fig. [5] No 



resonant effects are evident. Examination of the event distribution for m (tt + ti ) versus 
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Figure 5: Distribution of m(J/ipTT + ) for B° s — > J/ipir + iT~ candidate decays within ±20 
MeV of B° s mass shown with the (blue) solid line; m(J j '?fm + ) for wrong-sign J/?fm + 7T + 
combinations is shown with the (red) dashed line, as an estimate of the background. 



m 2 (J/?pir + ) in Fig. ^shows obvious structure in m 2 (ii + TT ) that we wish to understand. 



4.1.1 The signal function 

The signal function is taken to be the sum over resonant states that can decay into 7r + 7r~ 
plus a possible non-resonant S-wave contribution 



A=0,±1 



^afe^ l ^( Sl2 , S23 ,^; 



(3) 



where A^(si2, s 23 , 6j/^) is the amplitude of the decay via an intermediate resonance Ri 



with helicity A. Each Ri has an associated amplitude strength for each helicity state 
A and a phase The amplitudes are defined as 



A^s 12 ,s 23 ,e J/ ^) = f { b Lb) a r ( S23 ) f r Lr) tJ^-Y 

\TflB / 



p, 



R 



'S 2 3 



e 



(4) 



where Pg is the J/0 momentum in the B° s rest frame and Pr is the momentum of either 
of the two pions in the dipion rest frame, m# is the B Q S mass, F b Lb ' and F r Lr ^ are the B° s 
meson and Ri resonance decay form factors, Lb is the orbital angular momentum between 
the J /if) and 7r + 7r~ system, and Lr the orbital angular momentum in the 7r + 7r~ decay, 
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Figure 6: Distribution of S23 = m 2 (Tr + 7i ) versus Su = m 2 ( J/ipn + ) for B° s candidate 
decays within ±20 MeV of B° s mass. 



and thus is the same as the spin of the 7r + 7r~. Since the parent B° s has spin-0 and the 
J/ip is a vector, when the tt + tt~ system forms a spin-0 resonance, Lb = 1 and Lb = 0. 
For tt + 7c~ resonances with non-zero spin, Lb can be 0, 1 or 2 (1, 2 or 3) for Lb = 1(2) 
and so on. We take the lowest Lb as the default. 

<L R ) 



The Blatt-Weisskopf barrier factors F b Lb ^ and F A 
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arc 



F (0) 



1, 

y/z 2 + 3^ + 9 
\/z 2 + 3^ + 9 ' 



(5) 



For the B meson z = r 2 Pg, where r, the hadron scale, is taken as 5.0 GeV -1 ; for the R 
resonance z = r 2 P^ an d r is taken as 1.5 GeV -1 . In both cases z = r 2 P§ where Pq is 
the decay daughter momentum at the pole mass, different for the B° and the resonance 
decay. 

The angular term, T\, is obtained using the helicity formalism and is defined as 



(6) 



where d is the Wigner d-function |8j, J is the resonance spin, 9^ is the it + it resonance 
helicity angle which is defined as the angle of 7r + in the tt + 7t~ rest frame with respect to 
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the 7r + 7r direction in the B° s rest frame and calculated from the other variables as 

\m 2 (J/i)Ti + ) -m 2 (J/^TT-)]m(n + TT-) 

C ° S °™ = • (?) 

The J/ip helicity dependent term 0a(#j/^) is defined as 

®\{8j/i>) = J sin 2 Oj/^, for helicity = 
1 + cos 2 6 j/^ 

— ^ for helicity = ±1. (8) 

The function A R (s2s) describes the mass squared shape of the resonance R, that in 
most cases is a Breit-Wigner (BW) amplitude. Complications arise, however, when a new 
decay channel opens close to the resonant mass. The proximity of a second threshold 
distorts the line shape of the amplitude. This happens for the /o(980) because the K + K~ 
decay channel opens. Here we use a Flatte model [14]. For non-resonant processes, 
the amplitude v4 fl (s 23 ) is constant over the variables Si 2 and s 23 , and has an angular 
dependence due to the J /ip decay. 

The BW amplitude for a resonance decaying into two spin-0 particles, labeled as 2 
and 3, is 

M^) = — Ttt^ > ( 9 ) 

m z R - s 23 - im R T(s 2 3) 

where m R is the resonance mass, T(s 2 3) is its energy-dependent width that is parametrized 

as 

p \ 2L R +l , 



r(s ^H^J [vit) n - (10) 

Here To is the decay width when the invariant mass of the daughter combinations is equal 

tO TJlR. 

The Flatte model is parametrized as 

Ms™) = — : — r z ( n ) 

The constants g n7T and qkk are the /o(980) couplings to 7r + 7r~ and K + K~ final states 
respectively. The p factors are given by Lorentz-invariant phase space 




4ml 



m 2 (7i + Ti~ 



(12) 



il/'-^ + Wl-^- (13) 



m 2 (7T+7T ) 



The non-resonant amplitude is parametrized as 



P 

A(sta,S23,9j/^) = — J sin 2 Oji^. (14) 
rriB 
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4.2 Detection efficiency 



The detection efficiency is determined from a sample of one million B° s — > J/ip7i + 7i' 
Monte Carlo (MC) events that are generated flat in phase space with J/tj) — > ji + \i~ 



using Pythia 15 with a special LHCb parameter tune 16 , and the LHCb detector 



simulation based on Geant4 17 described in Ref 18 . After the final selections the 



MC has 78,470 signal events, reflecting an overall efficiency of 7.8%. The acceptance in 
cos 6 j is uniform. 

Next we describe the acceptance in terms of the mass squared variables. Both s\2 
and Sis range from 10.2 GeV 2 to 27.6 GeV 2 , where S13 is defined below, and thus are 
centered at 18.9 GeV 2 . We model the detection efficiency using the symmetric Dalitz 
plot observables 

x = s n - 18.9 GeV 2 , and y = s 13 - 18.9 GeV 2 . (15) 

These variables are related to S23 as 

S12 + si 3 + s 23 = m| + m 2 /v , + m 2 + + m 2 _ . (16) 

The detection efficiency is parametrized as a symmetric 4th order polynomial function 
given by 

£(s 12 ,s 2 3) = 1 + ei(x + y) + e 2 {x + yf + e 3 xy + e 4 (x + yf + e 5 xy(x + y) 

+e 6 (x + yf + e 7 xy(x + yf + e s x 2 y 2 , (17) 

where the Si are the fit parameters. 

The fitted polynomial function is shown in Fig. [7j The projections of the fit used 
to measure the efficiency parameters are shown in Fig. [8} The efficiency shapes are well 
described by the parametrization. 

To check the detection efficiency we compare our simulated J/ip<p events with our 
measured J/ip<p helicity distributions. The events are generated in the same manner 

1 1 2 

as for J/^7T + 7r~. Here we use the measured helicity amplitudes of L4||(0) = 0.231 and 
|A)(0)| 2 = 0.524 [5j. The background subtracted J/ip<p angular distributions, cos8j/^ and 
cos 9kk, defined in the same manner as for the J/ip7i + ir~ decay, are compared in Fig. [9] 
with the MC simulation. The x 2 / n df =389/400 is determined by binning the angular 
distributions in two dimensions. The p- value is 64.1%. The excellent agreement gives us 
confidence that the simulation accurately predicts the acceptance. 



4.3 Background composition 

The main background source is taken from the wrong-sign combinations within ±20 MeV 
of the B° s mass peak. In addition, an extra 4.5% contribution from combinatorial back- 
ground formed by J/ip and random p(770), which cannot be present in wrong-sign com- 
binations, is included using a MC sample. The level is determined by measuring the 
background yield as a function of 7t + tt~ mass. The background model is parametrized as 

B(s 12 , s 2 3, Qj/if) = Bi(s 12 , S23) x (l + a cos + /3 cos 2 6> J/V ,) , (18) 
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Figure 7: Parametrized detection efficiency as a function of S23 = m 2 {n + Ti ) versus 
S12 = m 2 (J / . The scale is arbitrary. 




Figure 8: Projections of invariant mass squared of (a) s\2 = m 2 (J/ip7i + ) and (b) s 2 3 = 
m 2 (7r + 7r _ ) of the MC Dalitz plot used to measure the efficiency parameters. The points 
represent the MC generated event distributions and the curves the polynomial fit. 



where the first part 5i(si 2 , S23) is modeled using the technique of multiquadric radial basis 
functions (19). These functions provide a useful way to parametrize multi-dimensional 
data giving sensible non-erratic behaviour and yet they follow significant variations in a 
smooth and faithful way. They are useful in this analysis in providing a modeling of the 
decay angular distributions in the resonance regions. Figure 10 shows the mass squared 
projections from the fit. The x 2 / n df of the fit is 182/145. We also used such functions 
with half the number of parameters and the changes were insignificant. The second part 
(l + a cos 0j/^ + (3 cos 2 0j/^) is a function of J / ip helicity angle. The cos 0j/^ distribution 

11 . fit with the function 1 + a cos 9 t/,i, + 8 cos 2 6 t/,i, that 



of background is shown in Fig. 
determines the parameters a = 



fit with the function 1 + a cos Oj/^p + 8 cos 9j/^ 
0.0050 ± 0.0201 and 8 = -0.2308 ± 0.0036. 
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5 Final state composition 



5.1 Resonance models 

To study the resonant structures of the decay — > J/tpn + n^ we use 13,424 candidates 
with invariant mass within ±20 MeV of the B° s mass peak. This includes both signal 
and background. Possible resonance candidates in the decay B° s — > J/^7r + 7r~ are listed 
in Table [TJ To understand what resonances are likely to contribute, it is important to 



Table 1: Possible resonance candidates in the B° s — > J/ip7i + 7i decay mode. 



Resonance 


Spin 


Helicity 


Resonance 
formalism 


/o(600) 








BW 


p(770) 


1 


0,±1 


BW 


/o(980) 








Flatte 


/ 2 (1270) 


2 


0,±1 


BW 


/o(1370) 








BW 


/o(1500) 








BW 



realize that the ss system in Fig. [T] is isoscalar (/ = 0) so when it produces a single meson 
it must have zero isospin, resulting in a symmetric isospin wavefunction for the two-pion 
system. Since the two-pions must be in an overall symmetric state, they must have even 
total angular momentum. In fact we only need to consider spin-0 and spin-2 particles as 
there are no known spin-4 particles in the kinematically accessible mass range below 1600 
MeV. The particles that could appear are spin-0 /o(600), spin-0 /o(980), spin-2 / 2 (1270), 
spin-0 /o(1370) and spin-0 /o(1500). Diagrams of higher order than the one shown in 
Fig. [I] could result in the production of isospin-one -k + it~ resonances, thus we use the 
p(770) as a test of the presence of these higher order processes. 
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Figure 9: Distributions of (a) cos 9j/^p, (b) cos9 KK for J/ip(f> background subtracted data 
(points) compared with the MC simulation (histogram). 
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Figure 10: Projections of invariant mass squared of (a) su = m 2 (J/ipir + ) and (b) S23 
of the background Dalitz plot. 



m 2 ('7r ± 7r ± 




Figure 11: The cos 9j/^ distribution of the background and the fitted function 1 + 

acosdj/^p + /3cos 2 Qjj^. 



We proceed by fitting with a single / (980), established from earlier measurements [T], 
and adding single resonant components until acceptable fits are found. Subsequently, we 
try the addition of other resonances. The models used are listed in Table [2] 

The masses and widths of the BW resonances are listed in Table |3l When used in 
the fit they are fixed to these values, except for the /o(1370), for which they are not well 
measured, and thus are allowed to vary using their quoted errors as constraints in the fits, 
taking the errors as being Gaussian. 

Besides the mass and width, the Flatte resonance shape has two additional parameters 
g-x-w and gxx, which are also allowed to vary in the fit. Parameters of the non-resonant 
amplitude are also allowed to vary. One magnitude and one phase in each helicity grouping 
have to be fixed, since the overall normalization is related to the signal yield, and only 
relative phases are physically meaningful. The normalization and phase of Jo (980) are 
fixed to 1 and respectively. The phase of /2(1270), with helicity = ±1 is also fixed to 
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Table 2: Models used in data fit. 



Name Components 



Single R /o(980) 

2R / (980) + / (1370) 

3R / (980) + / (1370) + / 2 (1270) 

3R+NR /o(980) + / (1370) + / 2 (1270) + non-resonant 

3R+NR + p(770) / (980) + / (1370) + / 2 (1270) + non-resonant +p(770) 

3R+NR +/ (1500) / (980) + / (1370) + / 2 (1270) + non-resonant +/„(1500) 

3R+NR +/ (600) / (980) + / (1370) + / 2 (1270) + non-resonant +/ (600) 



zero when it is included. All background and efficiency parameters are held static in the 
fit. 



Table 3: Breit-Wigner resonance parameters. 



Resonance 


Mass (MeV) 


Width (MeV) 


Source 


/o(600) 


513 ±32 


335 ± 67 


CLEO 




p(770) 


775.5 ±0.3 


149.1 ±0.8 


PDG 


rr 


/ 2 (1270) 


1275 ± 1 


185 ±3 


PDG 


|8| 


/o(1370) 


1434 ± 20 


172 ±33 


E791 


21 




/o(1500) 


1505 ± 6 


109±7 


PDG 


8 





To determine the complex amplitudes in a specific model, the data are fitted maxi- 
mizing the unbinned likelihood given as 

N 

£=n F fe s 23,^)> (19) 

1=1 

where N is the total number of events, and F is the total PDF defined in Eq. [T] The 
PDF is constructed from the signal fraction / sig , efficiency model e(si2, S23), background 
model B(s\2, S23, and the signal model S , (si 2 , S23, Ojm)- The PDF needs to be nor- 
malized. This is accomplished by first normalizing the J/ip helicity dependent part by 
analytical integration, and then for the mass dependent part using numerical integration 
over 500x500 bins. 



5.2 Fit results 



In order to compare the different models quantitatively an estimate of the goodness of fit 
is calculated from 3D partitions of the one angular and two mass-squared variables. We 
use the Poisson likelihood \ 2 22 defined as 

AW 



X 



2£ 



n, 



+ ndn 



71: 



X j 



(20) 
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where n.j is the number of events in the three dimensional bin i and X{ is the expected 
number of events in that bin according to the fitted likelihood function. A total of 
iVbin = 1356 bins are used to calculate the x 2 , using the variables m 2 (J/ip7i + ), m 2 {it + it~), 
and cos 6 j/^p. The x 2 /ndf and the negative of the logarithm of the likelihood, — ln£, of 
the fits are given in Table |4} There are two solutions of almost equal likelihood for the 
3R+NR model. Based on a detailed study of angular distributions (see Section 5.3) we 
choose one of these solutions and label it as "preferred". The other solution is called 
"alternate." We will use the differences between these to assign systematic uncertainties 
to the resonance fractions. The probability is improved noticeably adding components 



Table 4: % 2 /ndf and — ln£ of different resonance models. 



Resonance model 


-ln£ 


X 2 /ndf 


Probability (%) 


Single R 


59269 


1956/1352 





2R 


59001 


1498/1348 


0.25 


3R 


58973 


1455/1345 


1.88 


3R+NR (preferred) 


58945 


1415/1343 


8.41 


3R+NR (alternate) 


58946 


1414/1343 


8.70 


3R+NR + p(770) (preferred) 


58945 


1418/1341 


7.05 


3R+NR + p(770) (alternate) 


58944 


1416/1341 


7.57 


3R+NR + / (1500) (preferred) 


58943 


1416/1341 


7.57 


3R+NR + / (1500) (alternate) 


58941 


1407/1341 


10.26 


3R+NR + / (600) (preferred) 


58935 


1409/1341 


9.60 


3R+NR + / (600) (alternate) 


58937 


1412/1341 


8.69 



up to 3R+NR. Figure 12 shows the preferred model projections of m 2 (7r + 7r ) for the 



preferred model including only the 3R+NR components. The projections for the other 
considered models are indiscernible. The preferred model projections of m 2 ( J/-?/>7r + ) and 



cos9j/ip are shown in Fig. 13 for the preferred model 3R+NR fit. The projections of the 
other preferred model fits including the additional resonances are almost identical. 

While a complete description of the decay is given in terms of the fitted amplitudes and 
phases, knowledge of the contribution of each component can be summarized by defining 
a fit fraction, J-^. To determine J 7 ^ we integrate the squared amplitude of R over the 
Dalitz plot. The yield is then normalized by integrating the entire signal function over 
the same area. Specifically, 



j afe l ^^f(si2,s 2 3,6 l j/v) ds 12 ds 23 d cos 9 j/$ 



(21) 



/ S(si2,S23,9j/^) dsi2 ds 2 3 d COS 9j/tjj 

Note that the sum of the fit fractions is not necessarily unity due to the potential presence 
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Figure 12: Dalitz fit projections of m 2 (7r + 7r~) fit with 3R+NR for the preferred model. 
The points with error bars are data, the signal fit is shown with a (red) dashed line, the 
background with a (black) dotted line, and the (blue) solid line represents the total. The 
normalized residuals in each bin are shown below, defined as the difference between the 
data and the fit divided by the error on the data. 



of interference between two resonances. Interference term fractions are given by 

j.rr> _ ( J of a x e t{ ^~^ l4| : (si2, s 23 , d J/t )A^*(s 12 , s 23 , 9j/^)ds 12 ds 23 d cos Oj/,^ 
X \ J S(s 12 ,s 23 ,9 j/ip ) ds 12 ds 23 dcosOj/^, 

(22) 

and 

EfE^+E^'V 1 - ( 23 ) 

A \ R RR' / 

If the Dalitz plot has more destructive interference than constructive interference, the 
total fit fraction will be greater than one. Note that, interference between different spin- J 
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Figure 13: Dalitz fit projections of (a) su = m 2 (J/tpiT + ) and (b) cos #//,/, fit with the 
3R+NR preferred model. The points with error bars are data, the signal fit is shown with 
a (red) dashed line, the background with a (black) dotted line, and the (blue) solid line 
represents the total. 



states vanishes because the d J X0 angular functions in A x are orthogonal. 

The determination of the statistical errors of the fit fractions is difficult because they 
depend on the statistical errors of every fitted magnitude and phase. A toy Monte Carlo 
approach is used. We perform 500 toy experiments: each sample is generated according 
to the model PDF, input parameters are taken from the fit to the data. The correlations 
of fitted parameters are also taken into account. For each toy experiment the fit fractions 
are calculated. The distributions of the obtained fit fractions are described by Gaussian 
functions. The r.m.s. widths of the Gaussians are taken as the statistical errors on the 
corresponding parameters. The fit fractions are listed in Table [5] 

The 3R+NR fit describes the data well. For models adding more resonances, the 
additional components never have more than 3 standard deviation (a) significance, and 
the fit likelihoods are only slightly improved. In the 3R+NR solution all the components 
have more than 3<r significance, except for the ^2(1270) where we allow the helicity ±1 
components since the helicity component is significant. In all cases, we find the domi- 
nant contribution is S-wave which agrees with our previous less sophisticated analysis [3] . 
The D-wave contribution is small. The P-wave contribution is consistent with zero, as 
expected. The fit fractions from the alternate model are listed in Table [6} There are only 
small changes in the ^(1270) and p(770) components. 

The fit fractions of the interference terms for the preferred and alternate models are 
computed using Eq. 22 and listed in Table [7| 



5.3 Helicity distributions 

Only S and D waves contribute to the — > J/ipTT + ir~ final state in the m(7r + 7r~) region 
below 1550 MeV. Helicity information is already included in the signal model via Eqs. [7] 
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Table 5: Fit fractions (%) of contributing components for the preferred model. For P- 
and D-waves A represents the final state helicity. Here p refers to the p(770) meson. 



Components 


3R+NR 


3R±NR±p 


3R±NR±/ (1500) 


3R±NR±/ (600) 


/o(980) 


107.1 ±3.5 


104.8 ±3.9 


73.0 ±5.8 


115.2 ±5.3 


/o(1370) 


32.6 ±4.1 


32.3 ±3.7 


114 ±14 


34.5 ±4.0 


/o(1500) 


- 


- 


15.0 ±5.1 


- 


/o(600) 








4.7 ±2.5 


NR 


12.84 ±2.32 


12.2 ±2.2 


10.7 ±2.1 


23.7 ±3.6 


/ 2 (1270), A = 


0.76 ±0.25 


0.77 ±0.25 


1.07 ±0.37 


0.90 ±0.31 


/ 2 (1270), |A| = 1 


0.33 ±1.00 


0.26 ± 1.12 


1.02 ±0.83 


0.61 ±0.87 


p, A = 




0.66 ±0.53 






P, |A| = 1 




0.11 ±0.78 






Sum 


153.6 ±6.0 


151.1 ±6.0 


214.4 ±15.7 


179.6 ±8.0 


-ln£ 


58945 


58944 


58943 


58935 


X 2 /ndf 


1415/1343 


1418/1341 


1416/1341 


1409/1341 


Probability(%) 


8.41 


7.05 


7.57 


9.61 


Table 6: Fit fractions (%) of contributing components from different models for the 


alternate solution. 


For P- and D- 


waves A represents the final state helicity. Here p refers 


to the p(770) meson. 








Components 


3R±NR 


3R±NR±p 


3R±NR±/ (1500) 


3R±NR±/ (600) 


/o(980) 


100.8 ±2.9 


99.2 ±4.2 


96.9 ±3.8 


111 ± 15 


/o(1370) 


7.0 ±0.9 


6.9 ±0.9 


3.0 ±1.7 


8.0 ± 1.1 


/o(1500) 






4.7 ± 1.7 




/o(600) 








4.3 ±2.3 


NR 


13.8 ±2.3 


13.4 ±2.7 


13.4 ±2.4 


24.7 ±3.9 


/ 2 (1270), A = 


0.51 ±0.14 


0.52 ±0.14 


0.50 ±0.14 


0.51 ±0.14 


/ 2 (1270), |A| = 1 


0.24 ± 1.11 


0.19 ±1.38 


0.63 ±0.84 


0.48 ±0.89 


p, A = 




0.43 ±0.55 






P, |A| = 1 




0.14 ±0.78 






Sum 


122.4 ±4.0 


120.8 ±5.3 


119.2 ±5.2 


148.7 ± 15.5 


-ln£ 


58946 


58945 


58941 


58937 


X 2 /ndf 


1414/1343 


1416/1341 


1407/1341 


1412/1341 


Probability(%) 


8.70 


7.57 


10.26 


8.69 



and[8j For a spin-0 7r + 7r _ system cos9j/^ should be distributed as 1 — cos 2 9j/^ and cos 9^ 
should be flat. To test our fits we examine the cos9j/^ and cos9 niT distribution in different 
regions of 7r + 7r~ mass. The decay rate with respect to the cosine of the helicity angles is 
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Table 7: Fit fractions (%) of interference terms for both solutions of the 3R+NR model. 



Components 


Preferred 


Alternate 


/o(980) + / (1370) 


-36.6 ±4.6 


-5.4 ±2.3 


/o(980) + NR 


-16.1 ±2.7 


-23.6 ±2.6 


/o(1370) + NR 


0.8 ± 1.0 


6.6 ±0.8 


Sum 


-53.6 ±5.5 


-22.4 ±3.6 



given by 3 



dT 



d cos 9 j/^d cos 9 Tl 



A OQ ±^4 2O e^V5(3cos 2 0, 



sin 2 9j/^ 



(24) 



±- (|A 21 | 2 + 1 1 2 ) (15 sin 2 0™ cos 2 9J) (1 + cos 2 9 J/ij} ) 



where A 00 is the S-wave amplitude, A 2 i, i = —1,0, 1, the three D-wave amplitudes, and 
is the strong phase between A 00 and A 2 q amplitudes. Non-flat distributions in cos 9 n7T 
would indicate interference between the S-wave and D-wave amplitudes. 

To investigate the angular structure we then split the helicity distributions into three 
different 7r + 7r~ mass regions: one is the /o(980) region defined within ±90 MeV of the 
/o(980) mass and the others are defined within one full width of the / 2 (1270) and / (1370) 
masses, respectively (the width values are given in Table [3]). The cos9j/^p and cos0 T7r 
background-subtracted efficiency corrected distributions for these three different mass 
regions are presented in Figs. 14 and 15 The distributions are in good agreement with 
the 3R±NR preferred signal model. Furthermore, splitting into two bins, [—90, 0] and 
[0, 90] MeV, we see different shapes, because across the pole mass of / (980), the / (980)'s 
phase changes by it. Hence the relative phase between /o(980) and the small D-wave in 
the two regions changes very sharply. This feature is reproduced well by the "preferred" 

The "alternate" 



model and shown in Fig. 
description. 
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model gives an acceptable, but poorer 



5.4 Resonance parameters 

The fit results from the four-component best fit are listed in Table [8] for both the preferred 
and alternate solutions. The table summarizes the /o(980) mass, the Flatte resonances 
parameters g Tn , QkrIq-k-ki /o(1370) mass and width and the phases of the contributing 
resonances. 

The mass and resonance parameters depend strongly on the final state in which they 
are measured, and the form of the resonance fitting function. Thus we do not quote 
systematic errors on these values. The value found for the /o(980) mass in the Flatte 
function 939.9 ± 6.3 MeV is lower than most determinations, although the observed peak 
value is close to 980 MeV, the estimated PDG value (§]. This is due to the interference 
from other resonances. The BES collaboration using the same functional form found a 
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Figure 14: Background subtracted and acceptance corrected cosOj/^ helicity distributions 
fit with the preferred model: (a) in /o(980) mass region defined within ±90 MeV of 
980 MeV (x 2 /ndf =39/40), (b) in / 2 (1270) mass region defined within one full width of 
/ 2 (1270) mass (x 2 /ndf =25/40), (c) in / (1370) mass region defined within one full width 
of /2(1370) mass (x 2 /ndf = 24/40). The points with error bars are data and the solid 
blue lines show the fit from the 3R+NR model. 




Figure 15: Background subtracted and acceptance corrected cos 6*^ helicity distributions 
fit the preferred model: (a) in /o(980) mass region defined within ±90 MeV of 980 MeV 
(x 2 /ndf =38/40), (b) in ^(1270) mass region defined within one full width of ^(1270) 
mass (x 2 /ndf = 32/40), (c) in / (1370) mass region defined within one full width of 
/2(1370) mass (x 2 /ndf =37/40). The points with error bars are data and the solid blue 
lines show the fit from the 3R±NR model. 



mass value of 965±8±6 MeV in the J/ip — » <J)tt + tt~ final state [23]. They also found 
roughly similar values of the coupling constants as ours, g wn = 165 ± 10 ± 15 MeV, and 
Qkk/ 9-K-n = 4.21 ± 0.25 ± 0.21. The PDG provides only estimated values for the /o(1370) 
mass of 1200-1500 MeV and width 200-500 MeV, respectively [8]. Our result is within 
both of these ranges. 



5.5 Angular moments 

The angular moment distributions provide an additional way of visualizing the effects 
of different resonances and their interferences, similar to a partial wave analysis. This 



technique has been used in previous studies 24 



We define the angular moments (Y®) as the efficiency corrected and background sub- 
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Figure 16: Background subtracted and acceptance corrected cos^^ helicity distributions 
fit the preferred model: (a) in [-90,0] MeV of 980 MeV (x 2 /ndf =41/40), (b) in [0,90] 
MeV of 980 MeV (x 2 /ndf =31/40) 



Table 8: Fit results from the 3R+NR model for both the preferred and alternate solutions. 
indicates the phase with respect to the / (980). For the / 2 (1270), A represents the final 
state helicity. 

The parameters Preferred Alternate 



m fo (980) (MeV) 
ftnr(MeV) 

9kk/9 7T7T 

m fo (1370) (MeV) 
r /0 (i370)(MeV) 

0980 



; 1370 
^NR 
^1270, 
^1270, 



A = 
IAI 



939.9 ±6.3 
199 ± 30 
3.0 ±0.3 
1475.1 ±6.3 
113 ±11 
(fixed) 
241.5 ±6.3 
217.0 ±3.7 
165 ±15 
(fixed) 



939.2 ±6.5 
197 ±25 
3.1 ±0.2 
1474.4 ± 6.0 
108 ±11 
(fixed) 
181.7 ±8.4 
232.2 ±3.7 
118 ± 15 
(fixed) 



tracted 7r + 7r invariant mass distributions, weighted by spherical harmonic functions 

(Y®) = J dT(m nn , cos 6' 7r7r )l^ (cos O^dcos 9 nn . (25) 
The spherical harmonic functions satisfy 

/l r 
if(cos^)if (cos^tZcos^ = (26) 
.1 2tt 

If we assume that no 7r + 7r~ partial-waves of a higher order than D-wave contribute, 
then we can express the differential decay rate (dT) derived from Eq. (3) in terms of S-, 
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P-, and D- waves including helcity and ±1 components as 
dTim-xTTi cos 9 nn ) 

+ 2tt 



2tt |^ 5o F °(cos^) + ^e^oy^cos^) + Ao^Y^cose^ 

/ 2 

3 / 15 

A P±1 e t9p ^ \l — sm6 n7T + A D±1 e^ D ±^ J — sm8 n7T cos# 



8tt 



57T 



(27) 



where Ak x and <pk x are real- valued functions of m nn , and we have factored out the S-wave 
phase. We then calculate the angular moments 

A 2 So + A% + A 2 Do + A 2 P±1 + A 2 D±1 , 

4 [3 
2A So A Po cos 4> Po + —A Pq A Dq cos(0 Po - <j> Do ) + -A P±1 A D±1 cos(0 P±1 - 



4vr(F °) 

MY 3 °) 
4^(r 4 °) 



2v / 5 

-4 2 P „ + 2^ 5o> A Do cos0 Do + — A 2 Do - —A 2 



2 



1 



^ 2 



D±i> 



/ 3 6 

6 2 4 2 
jA Dq jA^ d±1 . 



Figure [17] shows the distributions of the angular moments for the preferred solution. 
In general the interpretation of these moments is that (Y®) is the efficiency corrected and 
background subtracted event distribution, (Yj°) the interference of the sum of S-wave and 
P-wave and P-wave and D-wave amplitudes, (Y" 2 °) the sum of the P-wave, D-wave and the 
interference of S-wave and D-wave amplitudes, (V 3 °) the interference between P-wave and 
D-wave, and (Y®) the D-wave. 

In our data the (Yj ) distribution is consistent with zero, confirming the absence of 
any P-wave. We do observe the effects of the ^(1270) in the (Y^ ) distribution including 
the interferences with the S-waves. The other moments are consistent with the absence 
of any structure, as expected. 



6 Results 
6.1 CP content 

The main result in this paper is that CP-odd final states dominate. The ^(1270) helicity 
±1 yield is (0.21 ± 0.65)%. As this represents a mixed CP state, the upper limit on the 
CP-even fraction due to this state is < 1.3% at 95% confidence level (CL). Adding the 
p(770) amplitude and repeating the fit shows that only an insignificant amount of p(770) 
can be tolerated; in fact, the isospin violating J/ipp(770) final state is limited to < 1.5% 
at 95% CL. The sum of / 2 (1270) helicity ±1 and p(770) is limited to < 2.3% at 95% CL. 
In the n + ir~ mass region within ±90 MeV of 980 MeV, this limit improves to < 0.6% at 
95% CL. 
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Figure 17: The tt + tc~ mass dependence of the spherical harmonic moments of cos 6 n7T after 
efficiency corrections and background subtraction: (a) (Yq), (b) (Y®), (c) (Y" 2 °), (d) (Kj°), 
(e) (Y®), (f) (Y®), (g) (Yq), and (h) (Y®). The points with error bars are the data points 
and the solid curves are derived from the 3R+NR preferred model. 
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6.2 Total branching fraction ratio 

To avoid the uncertainties associated with absolute branching fraction measurements, we 
quote branching fractions relative to the B° s — > J/ip<p channel. The detection efficiency 
for this channel from Monte Carlo simulation is (1.07 ± 0.01)%, where the error is due to 
the limited Monte Carlo sample size. 

The simulated detection efficiency for B° s — > J/ifm + 7r~ as a function of the m 2 (7r + 7r~) 



is shown in Fig. [L8j The simulation does not model the pion and kaon identification 
efficiencies with sufficient accuracy for our purposes. Therefore, we measure the kaon 
identification efficiency with respect to the Monte Carlo simulation. We use samples of 
D* + —7- 7T + D°, D° — > K~7i + events selected without kaon identification to measure the 
kaon and pion efficiencies with respect to the simulation, and an additional sample of 
K® — > 7r + 7r~ decay for pions. The identification efficiency is measured in bins of pp and 
t] and then the averages are weighted using the event distributions in the data. We find 
the correction to the J/ip4> efficiency is 0.970 (two kaons) and that to the J/ipfo efficiency 
is 0.973 (two pions). The additional correction due to particle identification then is 
0.997±0.010. In addition, we re-weight the B® p and p^ distributions in the simulation 
which lowers the tt + it~ efficiency by 1.01% with respect to the K + K~ efficiency. 

Dividing the number of the J/-?/>7r + 7r~ signal events by the J /%l)K + K~ yield, applying 
the additional corrections as described above, and taking into account B (0 — > K + K~) = 
(48.9 ±0.5)% |8), we find 

B (B° s J/^TT + TT-) , 

V ' — r-+ = 19.79 ± 0.47 ± 0.52)%. 

B (flo -> W) 

Whenever two uncertainties are quoted the first is statistical and the second systematic. 
The latter will be discussed later in Section [7j This branching fraction ratio has not been 
previously measured. 

6.3 Relative resonance yields 

Next we evaluate the relative yields for the 3R+NR fit to the J/ifj7r + 7r~ final state from 
the preferred solution. We normalize the individual fit fractions reported in Table [5] by 
the sum. These normalized fit fractions are listed in Table [9] along with the branching 
fraction relative to J/ip<f), <fi — > K + K~, defined as R r , where r refers to the particular 
final state under consideration. Thus 

B(B° s ^r) 

R r = ,- ■ 29 

B(B°^J/ij^ 

We use the difference between the preferred and alternate solutions found for the 3R+NR 
fit to assign a systematic uncertainty. Other systematic uncertainties are described in 
Section [71 

The value found for R r for the /o(980), 0.139 ± O.OOGI^JhiL i s consistent with the pre- 
diction of Ref. [9], and consistent with the our first observation using 33 pb -1 of integrated 
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Figure 18: Detection efficiency of — > J/tpir + 7r as a function of S23 = m 2 (-n- + 7i 



luminosity [IJ, after multiplying by B (<p — > K + K~). The decay B® — > J /ipj o(1370) is now 
established. Previously both LHCb jlj and Belle |2] had seen evidence for this final state. 
The normalized /2(1270) helicity zero rate is (0.49±0.16)% in the preferred model and 
(0.42±0.11)% for the alternate solution. 



Table 9: Normalized fit fractions (%) for alternate and preferred 3R+NR models and the 
ratio R (%) relative to B° s — >■ J/iptfi. The numbers for the /2(1270) refer only to the A = 
state. 



State 


Preferred 


Alternate 


R preferred 


R alternate 


Final R 


/o(980) 

/o(1370) 

NR 

/ 2 (1270) 


69.7 ±2.3 
21.2 ±2.7 
8.4 ± 1.5 
0.49 ±0.16 


82.4 ±2.3 
5.7 ±0.7 
11.3 ± 1.9 
0.42 ±0.11 


13.9 ±0.6 
4.19 ±0.53 
1.66 ±0.31 
0.098 ± 0.033 


16.3 ±0.6 
1.13 ±0.15 
2.23 ±0.39 
0.083 ± 0.022 


13.9 ±0.6±£§ 
4.19 ±0.53^ 
1.66 ±0.31+° ^ 
0.098 ± 0.033ig;g^ 



7 Systematic uncertainties 

Systematic uncertainties on the CP-odd fraction are negligible. In fact, using any of 
the alternate fits with different additional components does not introduce any significant 
fractions of CP-odd final states. 
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The systematic uncertainties on the branching fraction ratios have several contribu- 



tions listed in Table 10 Since R r is measured relative to J/ip(p there is no systematic 
uncertainty due to differences in the tracking performance between data and simulation. 
The J/ip<p P-wave yield is fully correlated with the S-wave yield whose uncertainty we 
estimate as 0.7% by changing the signal PDF, and the background shape. By far the 
largest uncertainty in every rate, except the total, is caused by our choice of the pre- 
ferred versus the alternate solutions. Using the difference between these fit results for the 
systematic uncertainty causes relatively large and asymmetric values. We also include 
systematic uncertainties due to the possible presence of the p(770), the /o(1500), or the 
/o(600) resonances by taking the maximum difference between the fit including one of 
these resonances and our preferred solution, if the difference is larger than the one be- 
tween the preferred and alternate 3R+NR fit. In the case of the /o(1500) the preferred 
solution is pathological in that it produces an unacceptably large / (1370) component 
along with a 214% component sum; therefore here we use the alternate solution that is 
much better behaved. 

The uncertainty from Monte Carlo sample size for the mass dependent 7r + 7r - efficien- 
cies are accounted for in the statistical errors, a residual systematic uncertainty is included 
that results from allowed changes in the shape due to the distribution of the events. The 
size of these differences depends on the mass range for the particular component multi- 
plied by the possible efficiency variation across this mass range. This is estimated as 1% 
for the entire mass range and is smaller for individual resonances. Small uncertainties are 
introduced if the simulation does not have the correct B° s kinematic distributions. We 
are relatively insensitive to any these differences in the B° s p and distributions since 
we are measuring relative rates. These distributions are varied by changing the weights 
in each bin by plus and minus the statistical error in that bin. We see at most a 0.5% 
change. There is a 2% systematic uncertainty assigned for the relative particle identi- 
fication efficiencies. These efficiencies have been corrected from those predicted in the 
simulation by using pion data from K® — > tt + tt~ decays and kaon and pion data from 

D*^ 1 —7- tt ± D°(D ), D°(D ) K^ir decays. The uncertainty on the corrections is 0.5% 
per track. The background modeling was changed by using a second-order polynomial 
shape in the J/^7r + 7r~ mass fit giving a 0.6% change in the signal yield. Since the input 
/o(1370) mass and width parameters were allowed to vary within Gaussian constraints, 
there is no additional uncertainty to account for. 

The effect on the fit fractions of changing the acceptance function is also evaluated. 
Since the acceptance model was tested by its agreement with the B° s — > J/ipK + K~ data 
in Fig. [9j we vary the data so that the model does not fit as well. This is accomplished by 
increasing the minimum IP x 2 requirement from 9 to 12.25 on both of the kaon candidates, 
which has the effect of increasing the x 2 / n df of the fit to angular distributions by 1 unit. 
The Monte Carlo simulation of B° s — > J/ilm + i{~ with the changed requirement is then 
fitted to get an acceptance function. This acceptance function is then applied to the data 
with the original minimum IP \ 2 cut of 9, and the likelihood fit is redone. The resulting 



fitted values from the preferred solution are compared with the original values in Table 11 
The changes are small and well within the statistical uncertainties. 
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Table 10: Relative systematic uncertainties on R(%). 



Parameter 


Total 


/o(980) 


/o(1370) 


NR 


f 2 (1270), A = 


m(n + 7i^) dependent effic. 


1.0 


0.2 


0.2 


1.0 


0.2 


PID efficiency 


2.0 


2.0 


2.0 


2.0 


2.0 


J/ip(f> S-wave 


0.7 


0.7 


0.7 


0.7 


0.7 


B° s p and px distributions 


0.5 


0.5 


0.5 


0.5 


0.5 


Acceptance function 





0.1 


1.3 


1.4 


3.9 


B(<j) ->■ K+K~) 


1.0 


1.0 


1.0 


1.0 


1.0 


Background 


0.6 


0.6 


0.6 


0.6 


0.6 


Resonance fit 




+18.2 
- 8.0 


+ 0.8 
-88.1 


+57.6 
- 3.7 


+ 3.0 
-15.8 


Total 


±2.7 


+18.3 
- 8.4 


+ 2.9 
-88.2 


+57.7 
- 4.8 


+ 5.5 
-16.4 



Table 11: Changes due to modified acceptance function. 



Values 


Original 


After change 


Variation(%) 


Fit fractions 








/o(980) 


(107.1±3.5)% 


107.2% 


0.1 


/ 2 (1270) A = 


(0.76±0.25)% 


0.79% 


3.9 


/ 2 (1270) |A| = 1 


(0.33±1.00)% 


0.26% 


21.2 


/o(1370) 


(32.6±4.1)% 


31.2% 


1.3 


NR 


(12.8±2.3)% 


12.7% 


1.4 


/o(980) parameters 








m fo (MeV) 


939.9±6.3 


938.4 


0.16 


ftnr(MeV) 


199±30 


205 


2.7 


9Kk/ 9tttt 


3.01±0.25 


3.05 


1.3 


/o(1370) parameters 






m fo (MeV) 


1475.1±6.3 


1476.4 


0.09 


T (MeV) 


112.7±11.1 


113.0 


0.27 



8 Conclusions 

We have studied the resonance structure of B° s — > J/ip7i + n^ using a modified Dalitz 
plot analysis where we also include the decay angle of the J/ip. The decay distributions 
are formed from a series of final states described by individual -k + -k~ interfering decay 
amplitudes. The largest component is the / (980) that is described by a Flatte function. 
The data are best described by adding Breit-Wigner amplitudes for the / (1370), the 
/ 2 (1270) resonances and a non-resonance contribution. Adding a p(770) into the fit does 
not improve the overall likelihood. Inclusion of / (600) or / (1500) does not result in 
significant signals for these resonances. 

Our three resonance plus non-resonance best fit is dominantly CP-odd S-wave over 
the entire signal region. We also have a D-wave component arising from the / 2 (1270) 
resonance. Part of this corresponds to the A 2 q amplitude which is also pure CP-odd and 
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is (0.49 ± 0.16io]o8)% °f the total rate. A mixed CP part corresponding to the A 2 ±i 
amplitude is (0.2 ± 0.7)% of the total. Adding this to the amount of allowed p(770), less 
than 1.5% at 95% CL, we find that the CP-odd fraction is greater than 0.977 at 95% 
CL. Thus, the entire mass range can be used to study CP violation with this almost pure 
CP-odd final state. 

The measured relative branching ratio is 

B (B° s ->• J/^TT + TT-) , 

v /_ — ^ = 19.79 ± 0.47 ± 0.52)%, 

B(B°^J/W) 

where the first uncertainty is statistical and the second systematic. The largest component 
is the / (980) resonance. We also determine 



B (B° s -> J/^tt+tt-) B (/ (980) -> 
B (1° ->• J lU) 



7T 1 7T 



;i3.9±0.6l^) 



This state was predicted to exist and have a branching fraction about 10% that of J/ip(j> [9]. 
Our new measurement is consistent with and somewhat larger than this prediction. Other 
models give somewhat higher rates 25]. We also have firmly established the existence of 
the J/'0/ o (137O) final state in P>° decay. 
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